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jxj '• 1 Introduction 

The neutrino emission processes play an important role in different astrophysical scenarios, 
in particular during the stellar core collapse and the cooling of a protoneutron star. Recently, 
several authors have studied some processes such as neutrino-electron scattering (Smit et al. 
1996| | and Cernohorsky, J. | |1994| |), neutrino coherent scattering off nuclei (Leinson [fL992[| ), 

or 



effects of nucleon spin fluctuations in the weak interaction rates (Janka et al. ||1996 



neutrino reactions with strange matter (Reddy and Prakash [1997]) in order to obtain the rates 



to be used in transport calculations. We will focus on the thermal pair emission-absorption 
process e + e~ «-> vv (TP in the next). In the stellar core collapse scenario, discrepancies in 
the efficiency of TP in the heating of matter outside the neutrinosphere have been noticed 
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(see, e.g., Janka, 1990 and references therein) and in order to obtain meaningful calculations 
of the energy deposition rate it is necessary to carefully consider the angular dependence of 
the distribution function. 

In some calculations involving neutrino transport, only the th and 1 st terms of the Leg- 
endre expansion of the collision kernel are included as in Bruenn ( ||1985|| ) and Suzuki ( ||1989|| ) . 



This approach could suffice where the diffusion approximation remains valid. The justification 
of that approach relies on the fact that the 2 nd and 3 rd order terms appear multiplying 2 nd and 
3 rd order terms of the Legendre expansion of the neutrino distribution function (I) and they 
vanish in the diffusion limit. However, if a general closure relation is used which is different 
to P = these higher order terms of / do not vanish, and their contribution becomes 
especially important in the semi-transparent region. In this paper we will also analyze the 
influence of the different closure relations that have been considered in recent years. 

This work is organized as follows: In §2 we present the Legendre expansion of the TP 
production and absorption kernels and we give explicit expressions for the th to 3 rd order 
terms. In §3 we study how these new terms affect the sources in the two-moment closure 
transport equations. In §4 we discuss the effects of the new terms and the influence of the 
closure relation. 

2 Legendre expansion of emission-absorption TP ker- 
nels 



Following Bruenn ([ |1985|| ) the contribution of thermal pair production and absorption of 
neutrino-antineutrino pairs to the right hand side of the Boltzmann equation (TP collision 
term) is 

B T p[I,I] = - . 3 / uj' 2 du' / dfi' / d(f) 
c(hc) Jo J-i Jo 

{[1 - 7][1 - I\I% P (w,u',coa9) -IIB% P (v,u',cosO)} (1) 

where I = I(t,r, fJL,u) (I = I(t,r, fi',u')) is the neutrino (antineutrino) invariant distribution 
function, u {uj') is the neutrino (antineutrino) energy in the frame comoving with the matter, 
H (//) the cosine of the angle between the neutrino (antineutrino) momentum and the polar 
axis, (ft is the azimuthal angle and 9 is the angle between neutrino and antineutrino directions. 
In what follows, the explicit dependence on t and r of the distribution functions will be omitted 
and we will assume axial symmetry with respect to the polar axis. The superscripts a and p 
refer to absorption and production, respectively. 

Assuming electrons and positrons in equilibrium at a temperature T, the following relation 
between the absorption and production kernels is satisfied 



R^p(uj,uj',cos9) = e^T R TP (u,u',cos9) (2) 
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where the temperature is measured in units of energy. Expanding R? TP as follows 

21 + 1 

B^ p (u,uj', cos 9) = $i(u,uj')Pi(cos9) 

i ^ 



(3) 



where Pi(cos8) are the Legendre polynomials and $>i(u,u') are the Legendre moments of the 
production kernel, the TP collision term can be written as 



B TP [I,I] 



2tt 



c(hcy Jo 



uj' 2 du' 



( 2/ -I- 1 



LO-\-U> 

1 — e T 



(4) 



In order to obtain the expressions for $/, we need to evaluate the following integral over 
electron energy E 



G 2 r+ J 



71 



r " dEF e (E, r])F e (uj + u/-E, -rj) x 
Jo 

ol\Ji(uj, lo', E) + ajJi(uj', u, E) 



(5) 



where ot\ and a 2 are summarized in table 1 for the different neutrino types and F e (E,rj) is 
the Fermi-Dirac distribution function 



UE,rj) = - s —- 
eT 'I _|_ 1 

rj being the electron degeneracy parameter, and 

1 



Ji(u,u',E) 



LOUJ 



7 9(u; W -E)x 



+1 



dfiPtin) A(fi) + B(ft)E + C(fi)E 2 9(/x-/i ) 



(6) 



(7) 



fiQ = 1 



2E(u + u)' — E) 



UOLO' 



In the previous expression O stands for the step function. 

Although A(n), B(fi) and C(/x) are complicated functions of /i, u and u', the integrals J/ 
can be done analytically, and the results for / = 0, 1 were first obtained by Bruenn ( ||1985|| ). 
We have calculated the integrals for I = 2, 3 and their dependence on E is simply a polynomial 
law of degree 2/ + 5 with coefficients being functions of u and to' . 
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Then, taking into account the following relation 



1 — e t 



F e (E, V )F e (uj + u J '-E,-r,) 

LO + a/" 



F e (E,r 1 )-F e [E,r 1 + 



T 



(9) 



the functions can be expressed in a simple way in terms of the dimensionless variables 
y = oj/T and z = u'/T. 



G 2 T 2 

% 1 - efe+ 2 ) 



[a 1 * l (y,z) + a 2 *i(z,y)] 



(10) 



^i(y, z) = Y^ (ci n G n (y, y + z) + d [n G n {z, y + 2)) 



n=0 



2i+5 



+ a in (G n (0, min(?/, z)) - G n (m&x(y, z), y + z)) 

n=3 



where 



G n (a,b) = dx ——— -j dx- 



X 



-(ri+y+z) _|_ I 



(11) 



(12) 



The explicit expressions for a; n , q n and <i; n coefficients are in Appendix A and the method 
to evaluate the G n (a, b) integrals is detailed in Appendix B. 

3 Energy and momentum transfer in two moment neu- 
trino transport 



Two moment neutrino transport (Cernohorsky and van Weert [|1992jD consists in solving the 
spectral energy and momentum balance equations as a coupled set. Let us define the I th 
moment of the distribution function 



1 r +1 1 



(13) 



and the following Eddington factors 



(14) 



Io(w) 



(15) 
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q(u) = 



/o(w) 



(16) 



For simplicity we will omit the energy dependence. Quantities with bar (7) will stand for 
antineutrinos and analogous definitions to (13-16) will be used. 

In order to close the set formed by the two equations (energy and momentum) we need 
two closure relations p = p(f,I ) and q = q(f,I ). In the next section we give some closure 
relations widely used in the literature. 

The source terms in the energy and momentum balance equations are obtained by angular 
integration of the collision term on the right hand side of the Boltzmann equation 



'dlo 
dt 



l r 1 



dfiB TP = 3 / uj' 2 duj' S (u,uj') 
-i c(hc) Jo 



TP 



2J-i 



c(hc) 



(17) 



'dh 
dt 



TP 



j d/ifiB T p 



l r l 

2 



2vr 



c(hc) 



/•CO 

/ J 2 dJ SAu,J) 
Jo 



(18) 



These source terms are the contribution of the TP process to the energy and momentum 
exchange between neutrinos and matter. After using the Legendre expansion (4) the expression 
for So and S\ can be obtained 



So 



'o - Jo 



e t 



(19) 



Sl = -l^ _ 7 1$1 + ( i _ e ^ ) Ti (lu, J) 



(20) 



r = E 



21 + 1$, f 1 



4 $ 



o J —1 J —1 



(21) 



(22) 



4 Discussion 

If neutrino transport calculations are done in the diffusion approximation, where the distri- 
bution function is truncated at first order 



I{n,u) = I (uj) + 3nh(uj) 



(23) 
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the expansion in (21) is also truncated at the same order and, hence, there is no need to 
calculate high order terms of the kernels. However, in the semitransparent region where the 
diffusion approximation breaks down and the flux factor is big, one needs to use a closure 
relation which is different from P = and the terms in <3>z for I > 2 must be taken into 
account because, as we will demonstrate, keeping only the first order term would give a wrong 
value of the energy exchange between neutrinos and matter. To see this fact more clearly, let 
us study the influence on the energy source of 2 nd and 3 rd order terms truncating r (c^, u/) at 
1 = 3. 



ToM = ! + {//} 3^ + { 



(3p-l)(3p- 1)^ , f( 5g-3/) (5g-3/) \ gg , . 

2 2 r$o 1 ' 




In this expression we have written inside brackets {} the part of the correction due to the 
distribution function. The value of these factors is restricted to the interval [0, 1]. Therefore, 
the maximum value of each new term is (21 + 1)J^, which we have plotted in figures la, lb 
and lc for I = 1,2, 3, respectively. As we can see from the plots, if terms in brackets are not 
small, the 2 nd contribution should be included for all energies and the 3 rd order term has a 
significant contribution for low energies. Only when the terms in brackets are much smaller 
than 1 the expansion can be truncated at first order. 

The closure relation becomes, therefore, a fundamental point in the calculation of the 
e + e~ <-» vv emission-absorption rate. There are several closure relations used by different 
authors and the question 'which is the best one?' has no answer yet. For the sake of com- 
parison, in this work we include four different closure relations: MB (Minerbo ||1978f |) , LP 



(Levermore and Pomraining ||1981| ), CB (Cernohorsky and Bludman ||1994| ) and MH (Mihalas 
|1984f ). Cernohorsky closure depends on both, the th and 1 st moments of the distribution 
function p = p(Io, f) and q = q(Io, f) and the rest of them are uniparametric closures, this 
is, p = p(f) and q = q(f). In the small occupation limit (Jo <C 1) CB closure is equivalent to 
Minerbo's one and in the maximal forward angular packing limit is equivalent to the vacuum 
approximation closure (see below). 

The form of these closures is the following: 

f(a) = coth(a) — 1/a 
MB { p(a) = l-2/(a)/a (25) 
q(a) = /(a) - (3p(a) - I) /a 



f(a) = coth(a) — 1/a 
LP { p(a) = coth(a)/(a) 

g(a) = coth(a)p(a) — l/3a 



(26) 



f p(/) = (l + 2/ 2 )/3 
U(/) = (3/ + 2/ 3 )/5 



(27) 
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CB {p(W) = § + f(l-/o)(l-2/ )x(i4) {2i 
I <i{Iqi f) — non-analytic 

where x{ x ) = 1 — 3x/j3(x) with (3{x) being the inverse of the Langevin function x = coth/3 • 
1/(3. The following expression fits this function well (Cernohorsky and Bludman [ 1994j| ) 



= = 



In figure 2 we show p(f) and g(/) for the different closure relations, taking an occupation 
level for the CB case of Iq = 0.1. In figure 3, the combinations of the different moments that 
appear inside brackets in (24) are plotted as a function of /. In both figures the solid line is 
for CB closure, dotted line for MH, dashed line for MB and dashed-dotted line for LP. We also 
plot with crosses the closure obtained in the vacuum approximation. As can be seen from the 
figures, there are important differences between different closures for / > 0.3. This can lead 
to relevant differences in the energy exchange between neutrinos and matter for large values 
of the flux factor /. 

To illustrate this feature, let us study a simple model consisting of a sphere of radius R 
radiating neutrinos and antineutrinos isotropically into vacuum, the so called vacuum approx- 



imation (Cooperstein et al. ||1986| ). The closure consistent with this model (VA) is 

f(a) = — ^ 

VA I p(a) = 1 + 3 + ° 2 (29) 

, . 1 + a + a 2 + a 3 
. = 4 

We assume that the neutrino (and antineutrino) spectrum at the surface of the sphere 
is Fermi-Dirac with zero chemical potential and T u = 1 MeV. Of course in a real case the 
neutrinospheres of neutrinos and antineutrinos are located in different places for different 
energies, but this example is illustrative of the general behaviour of the heating rate. We 
calculate the net (heating minus cooling) heating rate of the matter per unit volume for given 
distance to the center of the sphere (d) and matter temperatures (T). For T > T v cooling 
dominates over heating and we find that the effect of including new orders is negligible. In 
contrast, for low T the dominant term in eq. (19) is the term proportional to Tq and there are 
remarkable differences in the heating rates obtained including the new terms. These effects 
are also closure dependent and we compare the closure relations discussed previously in order 
to estimate the differences between them. 

In figure 4 we present the result of our calculations. We fix the matter temperature at 
T = 0.5MeV and we study the influence of higher order terms and closures for different 
distances from the center of the star (different flux factors). For the sake of comparison, 
we have performed Monte Carlo integration of the complete expression for the reaction rate, 
shown in the figure as the solid line. We overplot the total energy deposition after including 
each new order, using different symbols for the different closures. As can be seen, at first 
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order there is an underestimation of the deposition rate that is worse for high flux factors, 
even changing the sign (this means emission instead absorption of energy) for / > 0.8. The 
inclusion of the second order term gives a big improvement if one uses the closure consistent 
with the form of the distribution function in the model. The third order term is a small 
correction that can be omitted in all cases unless high accuracy is required. Using different 
closures instead of VA gives worse results at small x, but solves the problem in the sign for 
high values of x. We also observe remarkable differences between the results obtained using 
different closures, even though, we cannot deduce which one would have better behaviour in 
a realistic case. 

To summarize, we restate our main conclusions. The first conclusion is that, in the semi- 
transparent region and for matter temperature lower than neutrino temperature, it is necessary 
to consider the expansion up to 2 nd order. This procedure gives good results when combined 
with an appropriate closure relation. The 3 rd order term does not lead to a substantial 
improvement in the solution since, as we have shown, it is only a small correction. For T > T u 
cooling dominates over heating and we find that the effect of including new orders is negligible. 

The second conclusion is that, even though convergence is reached with 2 nd order correc- 
tions, the result is very sensitive to the closure relation chosen. Best results are obtained when 
using a closure relation consistent with the particular distribution function used in the model 
and therefore, detailed study of the closure in each particular problem is needed to obtain good 
estimates of interaction rates in a multigroup flux-limited diffusion problem. These results 
can be applied in all problems concerning neutrino transport such as stellar core collapse or 
cooling of newly born neutron stars. 
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Appendix A: Expressions of a/ n , c\ n and d\ n coefficients 

The coefficients for I = 0, 1 are the same as in Bruenn ( ||1985|| ) but adapted to our notation. 



Coefficients for / = 



Ay 2z 2 2y\ 
C00 = ^ ( — + yZ+ — > 

Az 3 

a, 



oo 



15y 



2 



coi = ~^( 3 ^ + 4z ) 
Az 2 

dm. - 



3y 2 

C02 



d, 



02 



a 03 



ao4 



a 05 



82 

3? 
8 

4 
3y 2 z 
A 



15y 2 z 2 
Coefficients for / = 1 

Ay .2 o 4 2 4 2 1 o. 
C10 = -- 1 ry 3 + ry 2 z + -yz 2 + -z 3 ) 
z 6 7 5 5 3 

4^ 3 



l io 



105?/ 



(142/ + 9«) 



%/ 4 2 7 2 2N 



Ay .3 1 

ci2 = -^ ( 5 y+ r J 



4z.l 3 , 

^12 = 3 + -z) 

j/" 5 3 5 



9 



^ - w 

4 

3y 3 z 



16 

a 17 



35y 3 z 3 
Coefficients for I = 2 



4^2 4 6 , 32 2 2 2 , 1 
-r(-Z/ + -2/ H 2/ ^ + -2/2 H 

4z 3 



105y 4 

4y.6 o 12 2 2 2 .. 
=-^(-y 3 + Y y 2 z + yz 2 + -z 3 ) 

42 2 48 , 6 1 



i21 = "^ ( 35 Z +7^+3^ 
8y 1 2 3 12 2 

8z ,1 , 3 12 2N 



023 " V 



4 

3y 3 z 



G24 = -^(10y + 9z) 



a 25 = 15^2 ( 7 % 2 + 126 ^ + 36 ^ 

12 . Q o 8 0\ 

a 2 6 = — 7-3(0 + 3yz+-z ) 
y 4 z 3 5 

48 

027 = — ^(2^ + 13^ + 12^) 



24 

a28 = ~W (2/ + 22;) 
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a29 " w 



Coefficients for I = 3 



4y 2 , 10 4 20 o 68 2 2 18 „ 3 4 , 

^30 = -^(V + 3% 2 ^ + 40yz* + 

4y 2 20 , 130 . 48 2 9 „. 
C31 = - ( 21^ + "63-^+35^ + 35^ 

4 * 2 /o 3 «, 2 130 2 200 3N 

dsi = 3V (% + % + ~ y + T" } 

V 50 2 6 6 2 
C32 = - — ( 63^ + 7 yZ+ 35 Z) 

Az 2 ,50 2 6 6 2 x 

cb = g- — z + + — y 

y b 63 7 35 



a33_ v 

4 , 

a34 = -T-3-(19y + 18z) 

a 35 = (253y 2 + 468yz + 180z 2 ) 

a 36 = -—^(U9y 3 + U7y 2 z + 330yz 2 + 50z 3 ) 
8 2Q1 fi 



ass = -l^E^y° + 54 ^ + 108yz 2 + 50/) 
40 , .„ 100 



a 3 (n) 



40 

9y 5 z b 
320 



a 3do) = -7^( 3 ^ + 52; ) 
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Appendix B: Evaluation of the G n integrals 

The integrals G n (a,b) appearing in the expression of the Legendre moments can easily be 
expressed as sums or differences of the Fermi-like integral 



F n (r),Xi 



dx 

n e x ~^ + 1 



and therefore, the problem is reduced to calculate this kind of integrals. 

In order to do this, we expand the denominator in the previous expression (Sack ||1990|| ). 
which must be done in a different way depending on the sign of a = x — rj, this is : 



1 



E(-i: 



a < 



m=0 



e a + l 



£(-l) m e-( m+1 ) a a>0 



m=0 



then, we obtain an infinite sum of integrals that can be calculated analytically using 



x k e mx 



1=0 



■l) k ~ l e mx x l 



I 



g 711 X 



1=0 



ttm k+1 - 1 



Let us define the following function 



Ti(a) = J2 



oo l j\n+l na 



n=l 



TV 



which is well defined for a > and / > 1, we finally arrive to a useful expression for the 
Fk(rj,xi) integrals depending on the value of rj. 



If T] < 



F fc (?7,xi) = k\ 



T k +i(-v) - 2^ — 



1=0 



l\ 



U0<rj<xx 



1 k+l 



k+l 



+ k\ 



INT\ 



1=0 



2 v T 2l+2 (0) V k -^ 
^ (k- l-2i)! 

T fc+ i_;(xi - Tj)x[ 



■l) k T k+1 (rj)-J2 



1=0 
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If X\ < T] 



k + 1 



k 



fc _;T fc+ i_;(7y ~ X^X' 
l\ 



1 1 

1 



E(-i) 



1=0 



The previous expressions are exact, and we only have to calculate a sum of a finite number 
of terms (up to k). The accuracy depends exclusively on the evaluation of the 7] (a) functions. 
The fact that these are uniparametric functions allows us to tabulate them in a fine grid at 
the beginning of the calculation and to obtain enough accuracy without excessive CPU time 
cost. 
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Figure 1: Ratios (2i + l)$j/$o for i — 1, 2, 3 as functions of y = ^ and z = ^ for an electron 
degeneracy parameter r\ e = 10 and \i-t neutrino type. 




Figure 2: Plots of p(f) and q(f) for different closure relations. The solid line is for CB 
closure, dotted line for MH, dashed line for MB and dashed-dotted line for LP. We also plot 
with crosses the closure obtained in the vacuum approximation. 
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Figure 3: Terms in brackets in the 2 and 3 order contributions in the Legendre expansion 
for different closures. The meaning of the lines is the same as in Figure 2. 
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Figure 4: Total energy deposited for different values of x — yl — (R/d) 2 in units of 10 20 erg 
cm -3 s _1 . The solid line is the exact solution after numerical integration and the different 
symbols stand for the closures LP (triangle), MH (star), MB (diamond) and CB (crosses). 



